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Abstract 

We directly compute Wilson loop/Wilson loop correlators on MxS 3 in 
AdS / CFT by constructing space-like minimal surfaces that connect two space- 
like circular contours on the boundary of global AdS that are separated by a 
space-like interval. We compare these minimal surfaces to the disconnected 
"double cap" solutions both to regulate the area, and show when the con- 
nected/disconnected solution is preferred. We find that for sufficiently large 
Wilson loops no transition occurs because the Wilson loops cannot be suffi- 
ciently separated on the sphere. This may be considered an effect similar to 
the Hawking-Page transition: the size of the sphere introduces a new scale 
into the problem, and so one can expect phase transitions to depend on this 
data. To construct the minimal area solutions, we employ a reduction a la 
Arutyunov-Russo-Tseytlin (used by them for spinning strings), and rely on 
the integrability of the reduced set of equations to write explicit results. 



1 Introduction 

AdS/CFT [1] has become an important tool for studying aspects of strongly coupled 
gauge theories (for a review, see [2]). In this context, classical minimal area world- 
sheets have played an interesting role: computing Wilson loop expectation values and 
quark/ant i- quark potentials [3-5]; computing the cusp anomalous dimension [6-8]; 
computing aspects of large spin or large R-charge operators [9-12]; describing the 
behavior of string vertex operator correlators for an AdS target space [13-16]; and 
computing color ordered scattering data [17] using the dual piecewise lightlike Wil- 
son loop. There has been great progress made in this area due to the integrability 
of the classical worldsheets on AdS (and sphere) backgrounds [11,12] as well as for 
the descendent Pohlmeyer reduced theories [18]. 

These computations often come down to evaluating a regularized worldsheet area 
and interpreting the exponentiation of this being some physical quantity of 

interest. We first define the Wilson loop operator as a path ordered exponentiation 
of the gauge field along a closed contour C 



The full M = 4 super Yang-Mills (SYM) theory Wilson loop is discussed in [4] and 
may contain some dependence on the scalar fields as well, however for our discussion 
we will will set this part to zero. The expectation value of the Wilson loop operator is 
identified with the appropriately normalized exponentiated classical area of a world 
sheet ending on the contours C at the boundary of AdS: 



Since we have set the scalar part of the Wilson loop operator to zero, we will be 
studying world sheets that have a profile in the AdS only. The solutions we will 
discuss are therefore applicable in any theory with an AdS factor in the geometry, 
however, one may need to be careful about the g s corrections to these backgrounds, 
as we shall discuss in a moment. 

In what follows, we will be concerned with two Wilson loop operators, defined by 
two circular contours. A natural object to study is [5,19], 



These have been previously studied for Minkowski space holographically in [19-22], 
and see [23] for a field theory discussion. One can argue for the above expansion 




(1) 



(W[C])=exp(-A[C}). 



(2) 




(3) 
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simply using large N arguments. The first two terms, 1 and -^2, are in some sense 
planar terms, where the 1 is enhanced due to being a "disconnected" diagram, and 
the order term comes from connected planar diagrams. Higher order terms come 
from non-planar diagrams. Further, we have extracted a certain power of A to make 
this look like a g s expansion, but is otherwise artificial, and could be absorbed into 
the definition of f n . 

In the string picture, we expect to compute a set of diagrams given in figure 1. 
The fact that we want to compute a ratio of the form (3) means that we want to 
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Figure 1: A diagramatic expansion of the worldsheets living in AdS. We have normal- 
ized the first term to be order one ((g s )°), so that the subsequent labels correspond 
to the degree to which they are sub-leading. 



regulate all of the areas of the above worldsheets by subtracting off the area of the 
disconnected leading order piece, i.e. to compute 

(W[Ci] 1 )\ij[C'l]) = 1 + 6XP [~ (^connected [@1 > ^] ~ ^disconnected [@1 > ) ] ( 4 ) 

where ^4c 0nnec ted * s the connected diagram with Euler character —In + 2. The dis- 
connected part is actually two distinct parts, coming from the area of the two dis- 
connected pieces. We recognize this expansion as being the same as the previous 
expansion because g s = gy M = 4, and the exponentiated area plays the role of the 
function f n . 

Here we have relied on the S 5 in the AdSsxS 5 geometry in a subtle way: we 
have assumed that the background has no g s corrections. Hence, for more general 
setups, like Sasaki-Einstein or orbifold backgrounds, more care may be needed; see 
[24] for recent work in this direction. In such cases, there is a second type of g s 
correction coming from the corrected worldsheet action due to the g s correction to 
the background. Therefore, lower genus diagrams will contribute to a given order 
in g s . For us, we will be discussing fx, and so this correction would amount to 
evaluating the old double cap solution in the correction to the sigma model action. 
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In fact, because we divide by a pair of (W) expectation values, we expect that this 
correction is in fact divided out and so we expect that the above is sufficient for 
leading order corrections in any background with an AdS. However, for higher 
corrections, one would need to be more careful, and check the contributions of lower 
genus worldsheets in the corrected action, as well as needing to modify the equations 
one solves. 

The purpose of this paper is to compute the first correction term in the above 
expansion in global AdS coordinates, namely, that associated with the "diagram" 
in figure 2. We compute this analytically for two synchronous coaxial same-size 

0X0-0) (D 

Figure 2: A diagramatic expansion of the first non trivial term we will compute. 
This gives the function /i at strong coupling. 

spacelike Wilson loops, separated by an interval in one of the angular coordinates 
((f)) of the S 3 inside of AdS 5 , and we find 

m (/l) = — (^connect — ^capx2) 

= -V^^l--^=[-(l-k 2 )K(k) + E(k)]^ (5) 
with k defined through the following relation 



/2ky/T^W ( \-k 2 „(2k 2 - 



sin 2 (# ) + tan 2 (^) + |tan (^) | cos(# ) 



sin 2 (fl ) + tan 2 (^) - |tan (^) | cos(0 o ) 



(6) 



and where 6q defines the "radius" of the two Wilson lines, and A0 determines the 
distance between their center points (see fig. 3). The functions K and II (and E 
later on) are complete elliptic integrals, see appendix B. The above is valid for k < k c 
where this critical value is given by 

- ^f=7 H 1 - k D K ( k c) + E(k c )YJ = (7) 
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Figure 3: Here we diagram the boundary conditions specifying the two equal size 
contours we use to explicitly calculate the Gross-Ooguri transition. The Cartesian 
limit gives the flat space behavior. 



which is numerically found to be 

k c = 0.8232. (8) 

Using this value of k c , one finds that the Gross-Ooguri transition happens at a A<p 
defined by 

A0 C \ sin 2 ^ 



tan 2 =P = " KU> (9) 

where £ ~ 2.4034 l . For A0 exceeding this value, the disconnected graph is preferred, 
and the transition is found to be first order, as in the flat Minkowski space case [19]. 
The flat Minkowski space limit is easy to read off, taking a. 9q limit, this reads 
Ax c = 0.9052i? agreeing with the result given in [19]. We also find that for sufficiently 
large Wilson loops that the connected solution always dominates, due to the finite 
separation distance possible on the sphere. This critical size is given approximately 
by 

sin 2 ^ rit ) = l {S ~ 1)2 



(£ + l) 2 

6 cvit ps 0.3547vr. (10) 

We see no problem with extending this to non- synchronous and non-equal-size Wil- 
son loops, given the structure of the equations we find in the subsequent sections, 



1 £ is just the value of cxp (^j=f= (^(fc) - ^j^-U (^i^, &))) evaluated at k = k c 
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although the above mentioned situation is the most symmetric. We comment on the 
structure of these more general surfaces in the text and in the appendix. 

One final note is in order. When the disconnected configuration is preferred, we 
expect that the correct correction is given by the disconnected area (i.e. 1, given 
that we have divided by this) times a propagation of a bulk field between the two 
disconnected caps: this is outside the scope of our current work, and has already 
been addressed in [5] for Minkowski space (for similar work for the null polygonal 
Wilson loops, see [25]). 

We organize the paper as follows. In section 2 we compute the equations of 
motion for the minimal surface in AdS n , and write the two dimensional (2D) Lax 
pair associated with this system. We reduce this Lax pair as in [11] and find a reduced 
one dimensional (ID) Lax pair. This furnishes a set of ID first order equations of 
motion, which we separate, and show how to solve by iterative integration. In section 
3 we write down several explicit solutions, including that for the synchronous coaxial 
Wilson loops, and synchronous coaxial equal-sized Wilson loops. We compute the 
regulated area minus the regulated area of the disconnected worldsheet to find the 
function j\ above, and compute the Gross-Ooguri transition point, and the critical 
radius where the transition no longer occurs. In section 4 we discuss possible future 
work. 

2 Equations of motion, reduction to ID system. 
2.1 Equations of motion, and reduction 

We first obtain the equations of motion for a string in AdS n . Recall that AdS n is 
given by the restriction 



in IR 2,n_1 . We will avoid using the flat space M 2 ' n_1 metric for computing the scalar 
product, and simply write We implement the above restriction via a Lagrange 
multiplier H in the Polyakov action, giving 



(y-i) 2 _ (Y ) 2 + (Y 1 ) 2 ■ ■ ■ + (y- 1 ) 2 



-1 



(11) 




(12) 



Upon elimination of S, the equations of motion read 



ddY 



M 



(dY ■ dY)Y 
-1 



M 



Y-Y 



(13) 
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and further, we have the Virasoro constraint 

dY ■ dY = By ■ By = o. (14) 

One may take the above (13) and antisymmetrize in a Y N to find 

Y [N d Q Y M] = Q ( 15 ) 

This is in fact invertible, one may obtain (13) from (15) by dotting in a Y , and 
using Y ■ Y — — 1. Finally, we rewrite the above as 

d (Y^BY M ^) + B (Y^dY^) = (16) 

to write the above as a total divergence. We can repackage these equations into 
equations of motion resulting from a Lax pair formulation, which we turn to now. 

First, writing the fields in some matrix g (which we explicitly show below), and 
defining A = dgg~ x ^A = Bgg~ l , the equations of motion can be obtained via the 
Lax pair with spectral parameter A 

a -lfl)* = - ( S -I^a)* = °- (17) 

The integrability condition for the two equations is 

( l-\)BA-(l + \)dA + \A,A] \ m 

(i+Aid-A) (18) 

which in the current context implies both the equation of motion (see footnote 2) 

BA + dA = (19) 

and gives the identity 

BA-8A+ [A, A] = 0. (20) 

In the above A, A plays the role of a gauge connection, and so we see that (17) is 
invariant under 

* = X®' 

A = xA'x-'+il + Xjdxx' 1 
A = xA'x^ + a-^dxx' 1 - 

(21) 
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Usually, we allow for any transformation xi, z t A) such that A! and A' are again 
independent of A. However here, our concerns will be more general and so below we 
will allow for any %■ 

Next, we wish to write these equations in Cartesian components, and for this, we 
recall that d a = |(9 + B),d T = \i{d — 8). Therefore, in these coordinates, we find 
the Lax pair 



da — 



.4 



.4 



2 Vl + A 1-A 
which can be written 



* = 0, 



d T — -i 



A 



A 



2 Vl+A 1-A 



* = 0. (22) 



dr-\ A \_ tX x ^ ))* = 0, (23) 



where A a = d a gg 1 ,A T = d T gg 1 . The gauge transformations can be read from the 
last ones quite easily and give 

$ = x^' 

K = xKx' 1 + d T xx~ l + ^d a xx~ l - 

(24) 

We define the combination ( Act 1 +^ 2 At ) = V a and similarly for V T . We will wish to 
reduce this model on a circle for a particular Ansatz, and we turn to this now. A 
particular choice of g for AdS 5 is given by 2 



(25) 



/ 







-z 3 


z 2 \ 




-Z\ 





z 2 


Zz 






-z 2 





-Z x 


\ 


-z 2 


-Zz 


Z x 


o / 



where we have defined 

Z x = Y- 1 + iY° 



Z 2 = Y 1 + %Y 2 = e i4> cos(9)r, 

Y 3 + iY 4 = e # sin(0)r. (26) 



2 Here we could have used the choice g = Y ■ T where T 1 are the SO(2,n — 1) (n=5) gamma 
matrices. This is particularly easy because dgg~ x = dYjYjT 1 J , and likewise for the 3 part of the 
connection. This makes the equation (19) read exactly as (16), so this is indeed the equations of 
motion. Here, to reduce the size of the matrices we are working with the 4x4 chiral components of 
the .50(2,4) matrices, which are a constant similarity transformation away from the g given here. 
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This gives a particular parameterization of the homogeneous coordinates Y 1 that 
solves the quadratic constraint (11). For this choice, the metric appears as 



ds 2 = riudYUY 3 

= -(l+r 2 )dt 2 + 



dr A 



1 + r 2 



(d# 2 + cos 2 (fl)# 2 + sin 2 (#)# 2 ) (27) 



We wish to make a reduction along a circle, and take that the a direction is 
compact, and we want to consider the solutions given by 



wa 



(28) 



and all other functions are functions of y only. Now, we ask the question of whether 
there exists some connection V a that is a independent after a gauge transformation. 
Indeed we find that there is one, and it is given by 



X 



(e 


-|tU(T 








\ 




























K 














(29) 



In fact, this makes both A a and A T independent of a altogether. One may be 
concerned that the functions above are not well defined on the circle (e.g. w an odd 
integer it is not well defined). However, going around the circle in a either takes \ to 
itself or \ to minus itself, but in the gauge transformation, \ only appears "squared" 
so that the resultant connection is well defined. The new connections A a and A T 
depend on A in some new way, however. Now we note that this simplifies our lives 
quite a bit. The Lax equations now read 



d T V a -d a V T + [V a ,V T } =0 
however, V T is a independent, and so we find 

d T V a = [V T ,V a ] 
which can be thought of as the one dimensional Lax pair 

U = [M, L) 



(30) 



(31) 



(32) 



where ' is now r differentiation. In such a formulation, the conserved quantities are 
given by Tr (L n ) because <9 T Tr (L n ) = nTr (L'L n_1 ) = nTr ([M, L]L n ~ l ) = by the 
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cyclicity of the trace. Since L is still a rational function of A, we expect there to be 
only a finite number of conserved quantities of this form. The above discussion is in 
fact the same reduction as was found in [11]: here we are using a subclass of Ansatz 
setting certain winding numbers to zero, and setting some dynamical functions to 
zero. However, here we will be considering the most general solutions to the system 
that we have found, rather than focusing on "constant radii" solutions discussed 
in [11]. 

In fact, the above Lax pair is insufficient to furnish all conserved quantities. 
However, it does furnish one that is not so obvious. To obtain a Lax pair furnishing 
all conserved quantities, one may simply redefine L + I(P t + \P<p) where I is the 
identity matrix, and P t and P^ are conserved quantities that we can find some other 
way 3 . This will not be needed for our further discussion, and we simply state a set 
of conserved quantities 

P t = (l + r 2 )t' 
P = r 2 cos 2 (0)0' 



H x = -(1 + r 2 ){t'f + + r 2 U6'f + cos 2 (0)(0') 2 - w 2 sin 2 (0)) (33) 

I _|_ r z 

H 2 = sin 2 (0)(l + r 2 )r 2 (t') 2 + ^Tg^Jg) + 2 s in(0) cos(0)r0V 

1 + H 

+ (cos 2 (0) - r 2 sin 2 (0)) r 2 {8') 2 + r 2 (- sin 2 (0) cos 2 (0)r V) 2 - w 2 sin 2 (0)) . 

Note that H\ is just the the hamiltonian resulting from the reduced action. Further, 
it must be for the above system to correspond to a minimal surface: Hi = is the 
Virasoro constraint. 

Next, we will eliminate t and from the discussion by replacing them with the 
conserved quantities P t and P^. Doing so, we find 

#i = -Txnr + T?~2 + r2 (( d ') 2 + ~ w2 sin2 (^)) ( 34 ) 

1 + r l 1 + r z \ r 4 cos / (6 l ) / 



P 2 r 2 sin 2 (0) (r') 2 sin 2 (0) 
H 2 = ^ + ^ ^ + 2 sin (6) cos V0V 



+ (cos 2 (0) - r 2 sin 2 (0)) r 2 (0') 2 + r 2 [ -3^^ - w * sin 2 (0) ) . (35) 

\ r z cos^(0J / 

The above quantities are conserved simply using the reduced equations of motion. 
However, since we are looking for a particular class of solutions, H\ — 0, we will find 



^Perhaps there is a more "efficient" Lax pair than the one we find here 
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other linear combinations more useful. We take instead the combinations 

(- cos(#)r' + sin(0)r(l + r 2 )9') 2 



(H 2 - E x 



1 + r 2 

[1 + r 2 sin 2 (#))((l + r 2 )P 2 - P 2 r 2 cos 2 (#)) 



n 2 



(36) 
(37) 



r 2 (l + r 2 ) cos 2 (9) 
(H 2 + r 2 sin 2 (0)#i) = (sin(0)r' + r cos(6)6') 2 
-w 2 r 2 sin 2 (6) (l + r 2 sin 2 (#)) 

where now the Virasoro constraint reads "Hi = % = "H- Note that 7^2 is only 
conserved quantities when the Virasoro constraint is met. Further note that the 
above combinations are special in that the single derivative term squared can be 
written as a total derivative. This allows us to write the above differential equations 

as 

r cos(6 l ) 



Or 



1 



\ 



■H + Pi + P? - Pi 



r cos(8) 

VT+7 2 



Pi 



l 



r cos(0) 
Vl+r 2 



(38) 

<9 r (rsin(fl)) = ± y/w 2 (r sin(#)) 4 + w 2 (r sin(fl)) 2 + H. (39) 

The second equation can be solve directly by integrating (note the right hand side is 
a function of rsin(0) only). Let us denote 

r cos(#) 



r sin(0), 



so that the above equations may be written 



F' = ±Jw 2 F* + w 2 F 2 + H 



l-F 2 



-n + p 2 + P 2 -p t 2 F 2 -p 2 T? . 



Along with the equations 



t' = P 



l-F 2 

i + f 2 



p. 



l-F 2 



F 2 (l + F 2 ) 
10 



(40) 

(41) 
(42) 

(43) 
(44) 



give all of the equations of motion. Therefore, one solves (41) for F2, then using this 
result, solves (42). This is possible because the right hand side of (42) is a function 
of F\ times a function of F2: since F2 is a known function of r, one can integrate 
this with respect to r and integrate with respect to Fi on the left, furnishing Fi 
as a function of r. With Fx and F 2 in hand, one may then directly integrates (43) 
and (44) find the general result. There is an important simplification here. When 
integrating the F\ equation, (42), one will be forced to do the integral J j^jv, and 

then F\ will be written as a function of this variable. However, this measure j^jv 
will appear when integrating to find t and <f> as well. Therefore, one can change to 
the coordinates defined by this integration measure, and simplify the integration. 
We will use this trick later, and emphasize it again when we turn to computing the 
solutions. 

Finally, if one wishes, one can solve the equations for non zero Hi, and this only 
changes the coefficient of _F 2 2 in (41), and the value of H appearing in (42), and so 
does not affect the method of solving the equations, nor the form of the equations. 
When considering profiles on the S 5 part of the geometry this will become important, 
but does not appear to be a major obstacle, given these simple modifications (see [26] 
for example). 

We now segue a bit and discuss the sections of F\ and F2 that are allowed. First, 
r is a positive definite number, and so the signs of Fi and F2 are actually given by 
9. Recall that for one cover of the S 3 inside of AdSs we have that < 9 < | to get 
a single covering. Extending 9 beyond this bound is relatively straightforward, as 
passing through 9 = can be compensated for by (— 9, (f>, if)) — > {9, (p, —if)), as these 
define the same homogeneous coordinates (and is the usual way for going through 
such an "origin"). Similarly, for going through 9 = | we see that (| + 9,(f>,if>) — > 
(I — 9, —(f), if)) . We can, therefore, make sense out of all signs for F\ and F 2 . However, 
note that — 1 < Fi = 77^=5 < 1- Hence, any solution that we find for F\ must lie in 
this range. Another way of saying this is that 

1 + F i A . Fil-Fl 



r 2 = ^2 - 1 , tan'(0) = ^f-£. (45) 



F{ V KJ F{ 1 + F< 



One should note that these expression define a real r and a real 9 for any real F\ , F 2 
that satisfy Ff < 1, and no restrictions on F2. 

There is one final note. Given the Virasoro constraint, one can get the pullback 
metric immediately as 

ds 2 = w 2 Fi (dr 2 + da 2 ) (46) 
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allowing direct computation of the lagrangian density as 

t 2 2 P 

S = t^L \ dTdoFl (47) 
4na' J 2 K J 

where we have reintroduced the AdS radius L. The above should be read with 
some care, however, when talking about Euclidean worldsheets: these require an 
extra factor of i. This will render our exponentiated action as being exponential 
suppression, and so should be understood as a tunneling event: we have already 
accounted for this in equation (2), where the area is to be read as the classical 
Euclidean area of the worldsheet. This exponential suppression is not surprising 
because we will be considering space-like separated 4 space-like Wilson loops. 



3 Solutions 

3.1 Special case: % = 

Let us work out a simpler set of equations as a warm up. We take % = above, and 
take ± = + for the F 1 equation and ± = — for the F 2 equation, and find 

F 2 (r) = —rr^ (48) 

sitin(w(r — T )) 



Fi(r) 



(P? ~ Pi) 2 (gg + (r - r ) - ltanh( W (r - r ))) 2 + g 
\ (Pi ~ Pi? (C 2 + (r- to) - itanhMr - r ))) 2 + P, 2 



(49) 



Note that this F 1 and F 2 define a good 9 and r given the restrictions mentioned 
above, and further, the signs were chosen so that < 9 < |. Particularly, F^ < 1 is 
always true as long as P t 2 > P?, i.e. that \P t \ > |P^|. 

Given the above functions, one can integrate the expressions for t(r) and </>(r) 
and further, one may rearrange the above expressions to solve for r(r) and 9{t). 



'All points of each Wilson loop are space-like separated from all points of the other Wilson loop 
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Upon doing so, we find 




cosh>(T - r )) ((P 2 - P 2 )<?o(r) 2 + ^ + ~? 



(w(r-r )) 



) 



(50) 



sinh 2 (w(r — r )) 



cosher - ro)) ((P, 2 - Pl?g a {rf + P|) 



(51) 



(52) 



(53) 



where we have defined the useful functions 



9o(i~) = C 2 + (r - r ) tanh(w(r - r )). 



(54) 



Note that in the above, we have the following constants of motions P t , P^, C 2 , to ; 00; r 0) 
and we have 4 dynamical fields. In general, we would expect 8 constants of motion 
for these equations, however, we must take H\ = to describe minimal surfaces, and 
we have taken the slice H = as well, hence we are left with six constants describing 
the surface. Finally, note that in the limit r — > that the surface goes to a circle, 
and t — > oo it goes to a point. This may have some bearing on the calculations 
performed in [27], however, we do not regulate these areas, and so we will have little 
to say on this point here. 

3.1.1 Sub-sub case: the single cap solution. 

We now explore a particular solution with Pi = const, % — 0, w = 1, Pt — P<j> = 
that will be important for later. Note that in the above discussion, plugging in 
"H = 0, one can take a limit with P 2 — > P^jP} being a fixed ratio less than one, but 
still get P^ = Pt = 0. In this limit, the solution reduces to 




T 



P!=P = COS(0 O ) 
t = t 

<p = 4> . 



(55) 
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where we have taken w — 1 for a single winding. For multiple windings, one can 
simply take r — > wt in the above equation. This "single cap" solution is well 
known [4,5] 

This gives a particularly simple solution 



COS 2 (fl ) + 7 



r 2 = f (T) (56) 

sin (6> ) 

9 = arctan ( *-^) . (57) 



This solution is the "single cap" solution. This is easiest to see in the following way. 
First, note that when r — > 0, the radius goes to infinity, and theta goes to some value 
fixed by 8 . Therefore, the constant sin(# ) determines the radius of the Wilson loop. 
Then, when r — > oo we see that 9 — > and r — > const, where again this constant is 
determined by 9 . However, we note that the pullback metric (46) gives 

ds 2 = \— (dr 2 + da 2 ) . (58) 

sinh (r) 

We would like to determine what the nature of the "cutoff" is when r — > oo. It is 
easy to see that in this limit, the pullback metric (46) reduces to 

ds 2 = e~ 2r (dr 2 + da 2 ) , (59) 

and changing variables to T = e~ T , we see that the metric becomes 

ds 2 = (dT 2 + T 2 da 2 ) , (60) 

where now the above limit is T — > 0. Because the periodicity of a is 27r, this represents 
a smooth point of the worldsheet, and is just a place where the round "cap" surface 
pinches off to zero size in a smooth way. If we had instead left w to be an integer 
winding number, the asymptotic form of the metric would have become 

ds 2 = w 2 e~ 2wT {dr 2 + da 2 ) (61) 

which can be transformed by the coordinate transformation T = e~ WT 

ds 2 = dT + w 2 T 2 da 2 . (62) 

This gives a conical singularity at the point T = 0, with an excess of angular "space" 
given by w. This, of course, is easy to visualize: the multiple windings are collapsing 
to a point, giving multiple covers of the flat space metric. 
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3.2 Finding F2 for general conserved quantities. 

Now let us look to the general H^O case. We write this as 



where 



(w 2 F 2 4 + w 2 F 2 +H)=H[1 




F 2 
1 2 



AH 

w z 



(63) 



(64) 



We recognize the above differential equation as being solve by Jacobi elliptic func- 
tions, however, it is not immediately clear how to write them in a manifestly real 
way. However, it is clear is to break the discussion into cases where certain reality 
conditions are met. We note that the combination y/1 — 7 is either real or imaginary 
depending on whether 7 is larger or smaller than 1 (w is assumed to be real). Hence, 
7 = 1 is a special value (in fact, this makes the right hand side of the equation 
above a perfect square, and this can be solved nicely too). However, there is also a 
crossover when — § + is positive or negative, in other words when 7 = (this 
is the case % = above). Hence we wish to break the problem into three cases: 

1. 7 < (H < 0) 

2. 0<7<1(0<^< w 2 /A) 

3. 7 > 1 (H > w 2 /A). 

After a bit of work, and using various identities [28-30], we can arrive at the 
following solutions for these cases 



±F 2 



(1-7) 



2 ' 2 



sn \ 2 
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7 ^ 

73 



7-1 



•dn I VVt' 72 7 3 



if 7 < 0, (H < 0) 

if < 7 < 1, (0 < H < 

if 7 >l, {%>*£) 



(65) 

where the relative "phases" associated with shifting r — > r — tq have been chosen so 
that the terms agree on the boundary of their relative regimes (7 = 0, 1). Above we 
have used a condensed notation where it is understood that all elliptic function have 
the same arguments, for example 

cn(z, k) 



cn 
sn ■ dn 



sn(z, k)dn(z, k) 



(66) 
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Figure 4: Here we have plotted F 2 (t) for one connected component of the worldsheet 
(where the arguments of the elliptic functions range between and 2K(k) for the 
appropriate k). The graphs depict different values of 7 = [—5,-2,-0.3,0,0.5,2,6] 
corresponding to the colors (line styles) [red (dot), orange (dash), brown (dash-dot), 
black (solid), dark green (wide space dot), blue (long dash), violet (wide spaced 
dash)]. 



Also note that in the above, the modulus (usually denoted "fc") of all functions is 
less than one. This gives some nice results. First, because < k < 1 we have that 
dn is positive definite, and both cn and sn are functions with amplitude 1. Thus, 
— is a generalization of a "-L" function. We can see this by remembering that all 

sn sin J 

elliptic function have period 4K(k) where K(k) is the complete elliptic function, and 
that sn(0, k) = and sn(2K(k), k) = 0: it goes to zero at and at the period over 
2. This then will have exactly the same shape as 4-. Similarly, — and are in 

J r sm J 1 sn sn-dn 

some sense generalizations of cot (blowing up at and period over 2, but crossing 
the horizontal r axis). This gives a rough idea of how these functions behave. We 
plot them to emphasize the point in figure 4. 

We note that in these graphs where % > that the worldsheet has been reflected 
through the origin: it has shrunk to zero size and then reemerged to finite size. We 
take that this is an orientation reversal of the world sheet, and this corresponds to 
reversing the orientation of the second Wilson loop operator. Hence, later when we 
define "coaxial" we will mean "coaxial with the same orientation." Reversing the 
orientation of one Wilson loop is considering the complex conjugate of the original 
Wilson loop operator, and we expect the correlator to depend on this information. 
We will not consider correlators of this form, "(H / [C]H^[C]''')" , and leave this to 
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future work. 



3.3 Synchronous coaxial Wilson loops, and the Gross Ooguri 
phase transition. 

We now turn to looking at a special class of solutions given a pair of Wilson loop 
operators that are synchronous, spacelike separated, and with the same orientation. 
Given these conditions and the above discussion, we see that this is a case where 
P t = 0, but Pff, 7^ and "H < (see also appendix A). In this case, the equations of 
motion simplify, and rather than arctan functions for the F\ integral, as in appendix 
A, we get instead logs. Solving these equations becomes relatively straightforward, 
as one simply needs to integrate to find Fx, and we find 



Fo(r) 



Fi(r) 



tr) 



7) 




7 



W + 4Pj + APiV, 



u 



(r) 



drP* 



1 -Fl 



W) 7 + 4Pi - APlH 



We will now turn our attention to the j^pi integral. For this case, we find 
dr 



(67) 

(68) 
(69) 




(70) 



IT ( sn(z- K(k),k), 
2k 2 

+n l,- 




where 



(1-7)" 



z = w(l — 7) 4 r 



(71) 
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and the definition of the elliptic pi function 5 is given in (130). The above answer 
is good enough for values of z in the interval (0,2K(k)), which is all we will really 
need (this interval defines a single connected component of the solution that begins 
and ends on the boundary of AdS). However, perhaps being a bit pedantic, we find a 
solution that is good for all values of z on the real axis (for the values of k needed). 
For this, we must be more careful with the signs of various functions, as II (T, a 2 , k) 
has branch-cuts when T = 1 for the values of a 2 and k given (both are less than one). 
For example, one should note that n(sn(z — K(k),k),a 2 ,k) behaves as a constant 
plus a square root of [1 — sn(z — K(k), k)} as sn(z — K(k), k) — > 1: note that sn(z, k) 
goes to 1 quadratically, and so the function is linear in z about this point. However, 
the sign of the linear correction is ambiguous: this is because we are evaluating y/z? 
which can be either ±x. Hence, we must put in a 0(sn(z, k)) to keep track of signs, 
where B is a step function that is +1 for positive values and —1 for negative values 
of its argument: 

{1 if x > 

undefined if x=0 (72) 
-1 ifz<0. 

However, now we see that we have introduced step function discontinuities at z — 
2mK(k) with integer m. These are easily removed with a "descending stairs" func- 
tion — % ln(e* x ) —x for some x given a branch cut for ln(z) that lies along the negative 

5 Here we have chosen to take a certain type of elliptic II function so that we avoid branch cuts 
at z = (2n + l)K(k), and instead have branch cuts at z = 2nK(k) to worry about. We do this 
because the later set of branch cuts occur at points where the worldsheet goes to the boundary of 
AdS, rather than in some interior point. 
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real axis (— oo,0). Putting this together, we find the final answer 



w 



dz 



(1-7) Wo l+^fc)) 2 ^ 

(1 - k 2 



(73) 
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fc e(sn(z,fc)) 
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K[k) 
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Ii-^n i, 



2k 2 — 1 \ 
n ( sn{z - K(k), k), - , fc ) e(sn(z, fc)) 



2 fc 2 -l 



• A- 



A- 



A; 2 ' 
In (e^^W" 1 ' 



P(A;) 



(74) 



One should keep in mind that the earlier function is sufficient for our needs, as this 
will describe one "boundary to boundary" section of the string. 
We have thus found the solution for P x given by 



Fxir) 



hir) 



C 2 exp 
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C 2 exp 
1-A: 2 ^ 



Jr/i ) + 4Pi - 4P 2 7* 
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2A- 2 - 1 

n ( sn(^ - K(k), k), — p— , fc ) e(sn(z, fc)) 



A- 2 



--In (V 7 ^" 1 ^- 



7T 



_L I 
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(1-7) 



1 > 
4 



z = w(l - 7) 4 r, C 2 > 



(75) 



and one takes the positive square root to determine Pi. We have included an "in- 
tegration constant" C 2 here, associated with the fact that we have defined p by a 
definite integral. This Pi lies in the acceptable range < P 2 < 1. This can be 
seen by noting that the numerator and denominator are both positive definite (given 
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that H < 0), and noting that the numerator minus the denominator is sign definite, 
given the sign of C 2 , and requires C 2 > so that it is negative. This restriction on 
C 2 is quite natural because it comes about as the exponentiation of an integration 
constant: if this constant is real, C 2 > 0. Finally, because F 2 is real, these functions 
define a real r and 9. 

The final integration for can be explicitly performed. This is because Fx is a 
function of the variable J j^pi- The integral actually has this factor j^ 2 ^ built in, 

and so we integrate this directly, obtaining 6 



— arctan 



( C 2 exp (^r/i) + 4P| + ±PpiS 



V / 



(76) 



So far, our discussion has been for general radius Wilson loops. In fact, one can 
see immediately that since Fx is in the appropriate range, and because Jj has no 
singularities at finite r, that Fi goes to constant values when F 2 diverges. This gives 

that the endpoints of < t < 2K< y k ) ^ are a ^ fa e boundary of AdS (this is the only 

10(1—7)? 

way to get a divergent F 2 and finite Fx). In fact, at these endpoints, C 2 helps fix the 
size of the two Wilson loops. 

Next, we note that the function I\ is monotonically increasing for increasing r: 
this is easy to understand because it is related to the integral J j^pi which has a 
positive definite integrand. Therefore we see that there is no way for the function 



g = C*2exp v 7 i Ix to be constant unless C 2 = 0, but this case is Fx — 1 which 

\(l-7) 3 / 

is forced to be on the boundary the entire time with = fixed. However, we may 
ask the question of when Fx has the same value at the two endpoints, namely when 

Fx(t = 0) = Fx ( r = K ^ ! I : this is where the two Wilson loops have the same 
V 10(1-7)*/ 



6 note that if we had let P t 7^ that the t equation of motion shares this property, so in the 
general case, both of these equations may be integrated using this variable to eliminate F% in a 
straightforward way. 
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size. This occurs when 



C 2 exp^/ I (0))+4^-4P|«y 

^ fe f ■ teo ) + 4 ^ 4 + 4 ^ 2w ) 2 - ^ (77) 
c2exp (^to)) +4i5 *- 4 ^ 2 



This happens when either 

; exp | ^r Ji(O) ) = Q>exp ( ( A(fc) ; ) 

(1-7)4 y ^(i_ 7 ) 4 ^(1-7)4 

which, again, due to Ji being monotonic, cannot happen unless C2 = 0; or they are 
equal when 



-0 



16P 4 (P|-H) : 



C 2 exp ( 7^-h-/i(0) = j VK 9 - j ' — ■ (79) 



(1—7)5 V w(l— T)i 
We have normalized so that I\(0) = and so we get 

Ci = 4P}(Pj-H) (80) 

exp ( -£5^h I -^r- 

V 2(1-7)5 Vw(l- 7 )3 



This value of C2 is special. It essentially dresses the factor of exp v \ I\{t) with 

\(l-7) 3 / 

a factor 4P| (P| — "H) and re-zeros the function Ii such that it is odd via reflection 
about t = K ^ 1 . This, in fact, renders the function Fi even with respect to this 

ui(l— 7) 5 

reflection. 

Therefore, the solution connecting two synchronous coaxial equal radius Wilson 
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loops is 
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(81) 



Due to the choice of C 2 one can readily check that Pi is symmetric with reflection 

about t = K ^ 1 , and further, so is P 2 . Finally, one can also see that the function 

10(1—7)3 



It, 



K(k) 



T,-P d 



\- const, and so the entire surface is symmetric 

210(1-7)* J 

under this interchange up to a choice of origin of 0, and so can be compensated for 
by shifting O - Note that the added reflection P^ — > — P^ does not affect the functions 
Pi or P 2 . Essentially all we are getting is that this very symmetric situation has a 
reflection symmetry that relates two halves of the surface (which is why we consider 
this case to begin with). Hence, if we wish to compare the connected surface to 
the disconnected surface, it suffices to do so for only 1 cap, and half the connected 
surface, which we turn to now. 
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3.3.1 Gross Ooguri transition. 



First, F\ at the boundary determines the radius of the Wilson loop. This is given by 

\ 2 

n 2 ( p |-^) 



Fx(r = 0) 2 = cos 2 (#(r = 0)) 



cxp 



(1-7)3 
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tu(l-7) 3 



+ 4Pi + 4P 2 ^ 
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cxp 



(l- 7 )3 



r 

0(1-7)3 



+ 4Pi-4P 2 H 



for the connected surface, and for the cap is given by 

Fl (r = 0) 2 = cos 2 (^ ) 



(82) 



(83) 



and so we have solved for the single-cap solution's one parameter, 6o, that gives the 
same boundary conditions (at this end of the connected surface). For the time being, 
we will concentrate on the case 

w = 1 (84) 

so that we are talking about a single winding, avoiding the conical singularity dis- 
cussed when constructing the cap solution. 7 In this case "H and 7 become redundant 
notation, and "H = 

Next, we need to define a regulator for both surfaces. Simply subtracting the 
answers is merely taking 00 — 00 and so can be adjusted to be any constant: the con- 
stant depending on the regulator. For this purpose, we introduce a radial regulator, 
cutting off both surfaces at r — -. For the single cap solution, this happens when 



t = arcsinh 




cos 2 (# )(l 




sin (0o J 



+ 



3 cos 2 



6 sin 



C(e E 



55) 



For the connected surface, we must consider the general expression 



F 2 



(86) 



7 However, the conical singularity may be OK in the current context. We are looking for the 
Gross Ooguri "phase transition." The argument is that once the minimal surface is actually the 
disconnected double cap type solution, the actual interaction is dominated by a bulk field. The 
site of the conical singularity may simply be an indication of where to put the vertex operator 
associated with this emission. 
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Note that if we expand the function for small values of r, we expect that F\ is 
approximately constant, and F 2 goes to 00 as -. Also, we note that I\ approaches 
zero as r 3 , and so we expect F± to approach it's constant asymptotic as r 3 . Also, 
the expression for Fi is an odd function, so we expect the series about r = to have 
an odd power series expansion. Hence, expanding the left hand side of the above 
expression as a series in r we should find 



% + C0 + O(r) = i (87 ) 



which we may solve order by order in e. For the first two terms, this is simple, and 
gives 

r = c_ 2 6 + ^pV + 0( e 4 ). (88) 

As argued above, these first two terms come from simply expanding F2 around r = 0, 
which we may do quite simply. Also, the asymptotic value of F\ matches the Fq = 
cos(# ) of the cap solution, so we find 

e 3cos^n)-i e3 + 0(e4) (gg) 



sin (#0 ) 6 sin 



3/ 



7 o> 



Indeed, this is because F 2 for the cap solution and F 2 for the connected solution 
agree for the first two orders in r when expanding around r = 0. Next, note that 
the Lagrangian density will diverge like -5, and so the action will diverge like -. 
The above prescription for the cutoff is accurate to a large enough order that the 
differences between regulating with the above cutoff and the radial cutoff will go to 
zero as the cutoff is removed. Hence, we use the cutoff 




r c = arcsinh — = (90) 

-cos 2 (# )(l + e 2 ) ' 

which is exact for the cap solution, and close enough for the connected solution in 
that deviations from the exact result will vanish as e — > 0. This argument suffices to 
show that we may do the integrals over the regulated range for r and then take the 
limit of the difference. 
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For this, we will need to have the integrals 
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(93) 



10 = 1 



where we have reintroduced the radius of AdS (L) and the appropriate factor of In 
from a integration, and E is the elliptic integral of the second kind. 
We may then take 



lim(yl conncct (e) 



A 



capx2 



(0) 
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[-(l-k 2 )K(k) + E(k)] 



(94) 



where E(k) is the complete elliptic integral of the second kind. This goes to zero 
at Ti — = 7. This actually serves as a bit of a check: recall that when % — > 
the actions coincide because the two different functions F 2 defining the two solutions 
coincide. 

There is in fact another zero, which we show by graphing the result in figure 5. 
Numerically, this value is given by y c ~ —6.92448860 which gives a critical value of 
k being k c = 0.82317491. 

Finally, to extract the other physical parameter defining the boundary conditions, 
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Figure 5: Above we have plotted the difference of the actions for the two worldsheets 
(red dashed line). We compute the critical value of 7 to be 7 C ~ —6.92448860 — > k c = 
0.82317491. We see that for values of 7 < 7 C that the connected solution is smaller, 
and so preferred, and for 7 C < 7 < that the disconnected surface is preferred. We 
will later translate this into a statement about the radius of the Wilson loops and 
the distance between them in <b. 



namely A0, we calculate 
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(95) 



which is rather odd looking at first, given that when P^ goes to zero, this quantity 
goes to infinity. However, this is just a statement about how to get cf> to go through 
a certain amount of angle. In some sense, at P^ = the string simply "drops" in 
the 8 or r direction without any rotation in 0. In this situation, the string passes 
very close to the place where dip 2 component of the metric goes to zero, and one 
picks up that <fi must go from to ir, i.e. that half of this change is n/2, which is 
indeed a tan(7r/2) = 00 limit. Hence, somewhat counter-intuitively, P^ = gives the 
maximum Acj) available. This, in fact, represents the maximally separated case for 
us because this is the only coordinate that separates the Wilson loops. What we are 
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finding is that if the two Wilson loops are separated by ir, the connected minimal 
surface passes very close to 9 = f , given that F 2 does not vanish at any point. 
Therefore, we have the two relations 
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(97) 



where on the left we have natural objects to define the boundary conditions F 
and A<j) and on the right, we have the two constants associated with the integrals of 
motion 7 and P^ (and recall that k is a function of 7). Now, we solve these equations 
and find 
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(100) 



where we have assumed a positive P^, which gives positive tan (4r). Note that the 
quantity on the right is strictly larger than zero, and larger than one. This serves as 
a check because I\ was defined to be positive definite for positive definite argument, 
and so the left hand side is strictly greater than one as well. To get a negative P^, one 
needs to simply reverse the sign of the right hand side of (98) because cos(# ) > 0, 
and further take the reciprocal of the right hand side of (99). Then, however, A</> 
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changes sign, so again, the new right hand side of (99) is bigger than one, again 
agreeing with the behavior of I\. This can be accounted for by simply taking the 
expression (99) and replacing tan (4r) — > [tan {-^) |, resulting in the expression (6) 
given in the introduction. In fact, the left hand side of (99) monotonically increasing 
for increasing 7 and the right hand side is monotonically increasing for increasing 
J tan (4r) | and fixed #0 (in fact, this expression is completely invertible). This gives 
that given < 8 < | that increasing 7 (decreasing it's magnitude, because 7 < 0) 
increases the A</> between the two Wilson loops. This gives the qualitatively right 
behavior. When the loops are close (7 — > —00) the connected surface is preferred, 
and when they are taken further apart (7 — > 0) the disconnected surface is preferred. 
This happens at 
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or solving 
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determines the critical separation distance where the disconnected surface is pre- 
ferred. One may take a flat space limit by considering small loops 9 — > 0. This 



gives 
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In this approximation, one can view the left hand side A0 as being the distance of 
a Cartesian coordinate, and the right hand side 9 as being the radius of the Wilson 
loops. Both sides have had units removed with the AdS radius. Hence, we interpret 
the above as saying 

(£-l) 2 
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£ 



-R z 



(104) 



where R is the radius of both Wilson loops, and Ax is the distance between their 
centers. Numerically, this gives Ax c = 0.90524981-R, agreeing with the flat space 
result quoted in [19]. 

There is one more interesting feature of the above equations. Given that £ is 
a constant, one may wonder whether there are sufficiently large Wilson loops such 
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that the critical distance A0 becomes maximal. Indeed, this happens when 

(l-^o,crit) = sin 2 (e cri t) = 1 ~ f £ ~ % ~ 0-8299619582 

9 CTit ps 0.36470575tt (105) 

and so these loops are actually quite large, becoming an appreciable fraction of the 
maximal angle 9 = We think of this effect as being similar to Hawking-Page [31], 
in that the scale associated with the sphere appears to play a crucial role. 



4 Discussion 

Here we discuss possible future directions. 

1. Note that we may describe "orientation reversed" Wilson loop correlators, given 
the structure of the solutions for 7i > 0. In these cases the string world- 
volume reflects through an origin, and reverses the orientation: this reverses 
the orientation of the second Wilson loop, and so computes a correlator of the 
form (WW''). This reverse orientation may have an interesting effect on which 
minimal surface is preferred, and so affect the Gross-Ooguri transition. 

2. It would be interesting to calculate the asymptotic behavior for small Wilson 
loops largely separated using the operator product expansion approach of [5] 
adapted to S 3 . This answer should agree with the flat space result [5] for a cer- 
tain regime, but should differ once the separation distance becomes comparable 
to the radius of the sphere. 

3. It would also be interesting to see how the transition [20] occurs for unequal 
size Wilson loops on S 3 , and see how this affects the "critical size" calculation 
done here for equal size Wilson loops. 

4. Having a profile on the S 5 should not change the behavior of the surfaces in 
AdS drastically: one can see this by noting the effect of treating H\ and H2 as 
both being non zero only affects certain coefficients in the equations of motion, 
rather than changing their structure. One expects the sphere equations of 
motion to completely decouple except for this effect (see e.g. [32] for a similar 
story for the gluon scattering minimal surfaces). 

5. One could also calculate correlators from insertion of other types of operators 
in the spirit of [5,33-37], and try to extract information about vertex operator 
correlators, as in [13-16]. 
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6. These solutions could also be used as a leading order "vaccuum" for soliton 
transformations. It would also be interesting to try to "add" multiple solu- 
tions together using generalized superposition principles sometimes available 
in integrable systems. 

7. It would also be interesting to see how finite size effects of the S 3 affect the 
behavior of small vs. large Wilson loops in the weak coupling approximation. 
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A Finding F\. 

A.l general considerations 

Now that we have F 2 , we may turn to finding Fx- First, we take an F 2 derivative of 
Fi by dividing the Fi equation by the F 2 equation to shed light on the final form of 
the solution: 



dF l 1 - Ff 



-H + Pl + P?-P?F?-P*-k 



dF 2 1 + F| yf% + W *F% + w 2 + Fi 
Therefore we need to perform two integrals 

dF 1 



(106) 



(1 - Ff)J-H + PI + Pt - PM - P 2 



dF 2 f dr 



Fi)^n + w 2 Fi + w 2 + Fi J (l + P 



h2\ 



(107) 



The F\ integral is what we would need to perform whether or not we had divided the 
equations, and the F 2 integral tells us the form of the final answer: it is some elliptic 
pi function (130). It is clear that the F 2 integral must be real when it is represented 
as a r integral, and so the integration for the Fx function must also yield a real result 
if Fx is to be real. 
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Therefore, let us check under what conditions one can find a real F\. It is quite 
clear that the Pi integral is only real if 

P f 2 + P%-H>0. (108) 

However, we can be more stringent simply by exploring the term under the square 
root more closely: 

-P?Ft + (P 2 + Pi - K)Fl - Pi (109) 

where we have multiplied by an P 2 for convenience (this has no effect on the sign of 
the above if Pi is real). We must have some sections where it is possible for Pi to 
be real: this happens where the above quantity is bigger than zero. This is easy to 
check simply by solving the above equation as a quadratic equation in P 2 . Solving 
this we find 



P 2 + PI-H- J (Pi + P 2 - UY - 4P*P* 



2P{ 



2 



<p'< 



p 2 + p 2 - n + A /(p t 2 + p 2 - uy - 4p 2 p 2 

(110) 



2P t 2 



p 2 +p$-h- y/m - p ) 2 - u) Wt + w - U) 



2P 2 " 1 " 



p 2 + p 2 - n + y/m - p<p) 2 - g) Wt + n) 2 - gj nin 

2P 2 ■ 1 ' 

One of these solutions must be real and positive for there to be a region where Pi 
can be real. In fact, because we know that P 2 + P| — H > all that needs to be 
shown (by the first line above) is that 

(P 2 + P 2 -H) 2 -4P|P 2 >0. (112) 

If this is the case, it is obvious that the + square root gives a positive value, but also 



the — value does also because y (P 2 + P| ~ %) 2 — 4P|P 4 2 is smaller in magnitude 
than (P 2 + P 2 — TL). The above inequality becomes 

\P? + Pl-H\ >2|P t ||P |, (113) 

however we have already seen that P 2 + P| — H > so the absolute value on the 
left hand side does nothing. This gives, finally, that 

(\P t \-\P^\) 2 >H. (114) 
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This restriction is more stringent than P t 2 + P? — Ti > 0, as seen from (113), and so 
the above restriction replaces it. The above restriction can, in fact, be read from the 
second line of (111) simply by considering the two cases where PtP^ is either positive 
or negative. In either case, one of the factors under the square root is positive, and 
so the sign under the square root is determined by the other factor: precisely the one 
that reads (\Pt\ — |P<^|) 2 — "H. We will see in a moment that these both are eliminated 
by an even more stringent constraint. 

There is one more important note to make: the function — 1 < Pi = r ^° s r ^ < 1. 
Above, we have just noted that for Pi to be real we must have Pf between the two 
solutions (111). The smaller of the two of these must be less than one to get an Pi 
that comes from a real r and real 9. Hence we find that 



p 2 + p$-h- J(p t 2 + pj - ny - ap»p? 

< i 



2P 2 



p*-h-J(p2-p?-h)*-4p?h 



-+ <0 (H5) 

We have already seen that the square root defined above must be real, via the 
constraint (|P t | - |P |) 2 > U. 

The above constraint breaks into two cases depending on the sign of %. When 
% > the first term dominates in magnitude compared to the term under the 
square root. This means that the first term determines the sign, so we find that 
P t 2 - Pi > —H. Combining this with (\P t \ - \P^\f > U we find 

|P t |-]P |>O for^>0. (116) 

This, therefore, tells us which square root to take from the previous inequality and 
we find 

\P t \-\P^\>Vn forH>0 (117) 

where above we take the + square root only. The TL < is trivial because the square 
root term in (115) dominates, and so the inequality is automatically satisfied. 

So, again, the above constraint is more stringent than all other constraints. Fur- 
ther note that the above is a strict inequality: if this number is equal to zero, then 
there is no window of allowed values for P^: this would give that P 2 = 1 strictly. 
This is only possible for r — > oo and so does not describe a true solution. In fact, we 
can see the inequality in the T-L — limit: 



h <a (118) 
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This has two distinct cases: if P? > P 2 then the above is equal to 0, and we get no 

solutions; if Pj < P 2 then we see that the above is t pi < which is consistent 

with the assumption. One can see this constraint from the exact solution given by the 
equation for 6 in (51). This also excludes the possibility that |P^| = \Pt\,l-i = 0: one 
can check this easily by plugging directly into the integrand finding f Fl =, 

which is manifestly imaginary. 

Finally, performing explicitly the Pi integral, one finds 

dF 1 



F?) X /-U + P$ + P?-P?F?-P^ 



arctanh f^S^HS^ 



Pl-P?+HHP?-Pj+K)Ff 



(119) 



up to a constant (that we leave in the P 2 integral). Consider first Ji < 0. Then, one 
can see that comparing the magnitudes of the numerator and denominator of the 
argument of arctanh 

{Pi -p? + n + (p 2 - p 2 + k)fI) 2 

-A(-H) (-P 2 P X 4 + (P 2 + P 2 - K)Fl - Pi) 

= {(P? + P} ~ nf - 4P?Pl) (1 - P 2 ) 2 > 0. (120) 

The result of the arctanh is therefore real because the argument lies in the interval 
(-1,1)- 

When Ti > 0, we see that we must push some factors of i through, which changes 
the function to an arctan. An arctan can accept any real input and so, other than a 
reality condition (which we have already met), the above is well defined. However, 
the output lies in a defined range ((— | , |J, and one may think that this constrains 
the value of the integration constant in P 2 . We will address this when the time 
comes, and find that it is not an issue. 

There is an extremely important caveat for the above arguments to hold: that 
we have not divided by a function that is identically 0. This case deserves special 
consideration. We can immediately see that if this term is zero, we have a constant 
Pi. In such a case, we find that the 



p 2 + P!-n± J(p 2 + p 2 - ny - 4p 2 p 2 



F't = — : " -o — : : • (121) 

1 2P 2 v ; 



33 



For this, we have that one of these solutions must satisfy F 2 < 1 as well, so that the 
above restrictions found must also hold. However, we must also explore the + square 
root, and see when it is possible for this to be less than one. We therefore explore 



P 2 + PI--H+ J(P? + P 2 - Uy - 4P 2 P 2 
°< y ~2Pf 



p 2 + pI-h + a /(-p 2 + p 2 - uy - AP?U 



-1 < ^2 < 0- (122) 



2J? 



When H < 0, the bottom line is clearly violated on the upper bound (given a real 
square root). When H > 0, the top lines lower bound is obviously satisfied, and the 
lower lines upper bound is satisfied if P 2 — P? + % > 0: this condition is the same 



as when the square root is the other sign, and so \P t \ — |P^| > sfH. 
To conclude, we summarize the restrictions found in this section 

For F[ ^ 0: % < : No restrictions 

n>o-. \p t \- ipj > Vn 



For F[ = 0: P x 2 



P 2 + P 2 - U ± ,/(P t 2 + P, 2 - ?/) 2 - 4P 2 P 2 



2P 



H < : minus sign only, no other restrictions 

ft > : \P t \ - |PJ > y/H either sign OK (123) 

Note that the other possible solution P 2 = r = 1 for P{ = gives that the 

string is exactly on the boundary, and so this is not an interesting solution. 

A.2 Finding F 1 for U < 0, 7 < 

Next, let us name a function g(r), given in 119 



2V=HJ-P?F} + (P 2 + P 2 - U)Fl - PI 



P*-If + H + (I?-P} + K)F* 

dr 

'2 



g(r) = - tanh ( 2^H \ 7^ ) (124) 
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where the P2 integral is solved the same way as in the text. One may solve the above 
for F\ to find 



m - p lf - K 2 ) 9(r) 2 - 2H(P 2 + Pi - U) 



±2U^{1 - g{rf) ((P 2 + P 2 - - AP 2 Pf) 
+ [(P?-P* + n) 2 g(T) 2 -4P t 2 U} 



(125) 



One may check that the expression in strictly greater than zero because the term 
under the square root has smaller magnitude than the term without a square root 
(the difference of squares is positive) for the allowed values of < g(r) 2 < 1. Next, 
one may compute 



F 2 - 1 



- 2n(P 2 - Pi + n)g(r) 2 - 2H{Pl - P 2 - H) 



±2UsJ{i - 9 {tY) ((p 2 + p| - uy - ap?pd 

^[{P 2 -Pl + V) 2 g{r) 2 -AP 2 U\. 



(126) 



In this, one can show that the term with the square root dominates, and so this term 
completely determines the sign. Therefore, we must take the + sign above (because 
U < 0) to have F? < 1. Hence, we find 

F 2 = ((P 2 - Pi) 2 - U 2 ) g{r) 2 - 2H{P 2 + P 2 - H) 

+27^(1 - g{rf) ((P t 2 + PI - Uf - ApJpJ) 

- [(P 2 - Pi + H) 2 g(r) 2 - AP 2 U] (127) 
where this expression is manifestly positive and less than 1 for < g(r) 2 < 1. 
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B Elliptic Functions 



We define the elliptic integrals 



r l 

F(T, k) = / dt , 



(128) 



E(T, k) = [ dt 
Jo 




(129) 




o (l-a 2 t 2 )y/(l-t 2 )(l-kH 2 ) 



(130) 



and the complete elliptic integrals as 



K(k) = F(l,k) 
E(k) = E(l,k) 
U(a 2 ,k) = 11(1, a 2 , k). 



(131) 
(132) 
(133) 



For the Jacobi elliptic functions, we follow the conventions of [28-30]. The refer- 
ence [28] has several forms of the above integrals useful for defining various com- 
binations of the Jacobi elliptic functions, while [29] has extremely useful tables for 
their transformation properties under shifts, and [30] has some useful transformation 
properties of II (T, a 2 , k) not contained in the other references. 
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